Condensation of Ideal Bose Gas Confined in a Box Within a Canonical Ensemble 
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We set up recursion relations for the partition function and the ground-state occupancy for a fixed 
number of non-interacting bosons confined in a square box potential and determine the temperature 



dependence of the specific heat and the particle number in the ground state. A proper semiclassical 
treatment is set up which yields the correct small- T-behavior in contrast to an earlier theory in 
Feynman's textbook on Statistical Mechanics, in which the special role of the ground state was 
ignored. The results are compared with an exact quantum mechanical treatment. Furthermore, we 
£Nj i derive the finite-size effect of the system. 



PACS numbers: 03.75.Hh 



I. INTRODUCTION 



The thermodynamic properties of a Bose gas are usually studied in the grand-canonical formulation where energy 
and particle number are fixed on the average. Present-day experiments, on the other hand, are performed in a 
magnetic or optical trap which contain a roughly fixed number of particles. Thus, the experimental situation may 
be better approximated by a canonical ensemble with fixed particle number. In the thermodynamic limit one usually 
expects grand-canonical and canonical treatments to yield the same results for all thermodynamic quantities and 
the condensate fraction. This has indeed been proved in [l[ in general and for harmonic traps in 0, y, 0, 0). Some 



statistical quantities, however, do not have the same thermodynamic limits in the two ensembles. One fundamental 
example is the fluctuation width ANq = W ((ANq) 2 ) of the particle-number No in the ground state. In the grand- 
canonical ensemble, this quantity has the size \J Nq{Nq + 1) and is at low temperatures T of the order of the total 
' particle number A" 0, 0, M S S]- In the canonical ensemble, on the other hand, this quantity can be shown to 
vanish at T — since No becomes equal to the fixed particle number N. The grand-canonical result is a consequence 
of the assumed statistical independence of all open grand-canonical subsystems. This assumption is only fulfilled if 
£ — . there is no long-range ordered condensate, otherwise it is wrong. The problem can be properly removed by fixing the 
whole particle number and ignoring its fluctuations as was shown e.g. in [H, 0]. The particle number in the ground 
f*"*. ■ state is then Ao = N — N cx , where AT CX is the excited particle number. Since, furthermore, the fluctuation AN CK 
of the thermal gas, consisting of excited states, can be shown to vanish for zero temperature, the same is true for 
AAo. Its temperature behavior is found to be AAo ~ (T/T c ) 3 / 2 A 1 / 2 for an ideal gas in a harmonic trap 0] and 
AN ~ (T/T c ) N 2 ^ for a homo geneous system [9( for temperatures T below its critical value T c . By fixing the full 
particle number AT, the system does not exchange its particles with the environment at all and this situation, therefore, 
corresponds to one in a canonical ensemble. The only exception is that excited states are then still described by the 
grand-canonical Bose-Einstein distribution. The entirely canonical description will be presented below. 

The canonical treatment is conveniently performed in a first-quantized path-integral description. All particles are 
represented by cycles winding around a cylinder in Euclidean spacetime whose circumference is the imaginary-time 
axis t € (0, h/ksT) 13, 11 1- For sufficiently low temperatures, the indistinguishability of particles leads to long closed 



paths winding many times around this cylinder, giving rise to correlated subsets of particles [12]. In Ref. [131 ]. this 



representation was utilized to claim that condensation amounts directly to a proliferation of long cycles. This claim 
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can, however, not be upheld. The formation of correlated collectives is the origin of superfluidity, not condensation 



as it was correctly implemented for an ideal gas in a harmonic trap [14[. The essential signal for the Bose-Einstein 
condensation is, instead, associated with the occupancy of the single-particle ground state. Thus, not the winding 
numbers of the cycles are relevant to condensation, but the weight of the ground state in them. For a harmonic trap, 
this was properly shown in Refs. ljl 16 1. 



If the harmonic trap is replaced by a square box, whose infinite-volume limit is the homogeneous system, the 
systems confined in it show some unusual features. It is the purpose of this paper to exhibit these. The discussion 
starts in Section HT1 with a brief summary of the imaginary-time evolution amplitude for a fixed number N of bosons, 
where the indistinguishability of particles is explicitly taken into account. In Section IlIII we calculate the partition 
function of such a canonical system from the trace of the TV-particle imaginary-time evolution amplitude written in 
the cycle representation of paths. Since this is hard to evaluate numeric ally for a large particle number, we have 
derived an efficient recursion relation along the lines of Refs. [1, [l?], 18, 19, H3] in Section ITVl In Section fVl this 
relation is made specific by a spectral decomposition of the canonical partition function. This allows us in Section IVT1 
to extract the ground-state contribution to the partition function from the cycle representation, from which we read 
off the probability to find a particle in the ground state. This is used to define the condensate fraction Nq/N for a 
system with a fixed particle number TV. In Section TVI!! we apply our results to the Bose gas confined in the square-box 
potential. There we show that the semiclassical approximation requires a delicate treatment of the ground state for 
finite systems. This is necessary in order to obtain the correct low-temperature behavior of the thermodynamical 
functions in contrast to an earlier theory in Feynman's original textbook [lfj. The results are then compared with 
an exact quantum mechanical calculation. In the subsequent Section IVIIII we investigate the finite-size effect in a 
box potential. We identify a would-be transition temperature as the finite- N version of the critical temperature of 
Bose-Einstein condensation in a canonical ensemble. Subsequently, canonical results calculated for these crossover 
temperatures arc compared with grand-canonical ones. We end with an analytic calculation of the finite-size effect 
upon the critical temperature in the grand-canonical ensemble. 



II. QUANTUM STATISTICS OF IDENTICAL PARTICLES 



We start by briefly reviewin g th e quantum statistics of a fixed number N of non-interacting identical particles in 
a first quantized approach [13, U|- At first, they are treated as distinguishable particles. Their orbits in imaginary 
time r are denoted by x„(r) with the particle indices v = 1, 2, . . . , N. Global and local thermodynamic properties 
of the many-particle system are determined by the imaginary-time evolution amplitude which is given by the path 
integral 



N 



(x 



16) 



> x JVb, Tfe|x la , . . . , Xjv a , T a ) = < Jj[ 



x„(t,,)=x„ 



D X „ \ e --4[xi,...,x^]/fi _ 



_„ = 1 •/x„(r a )=x„ 

The Euclidean action in the exponent has the generic form 



-A[xi,. 



N 

E 



M 



±l(r) + V(Mr)) 



(1) 



(2) 



with the particle mass M and the common background potential V (x„) for all particles. Due to the additivity of 
the action ((2|) we find that the imaginary-time evolution amplitude of N distinguishable non-interacting particles 
factorizes into N one-particle amplitudes 



(xifc, . . . , XjVb, Tb|xi Q , . . . , XjVq, T a ) = (Xi b , T b |xi Q , T a ) ■ • ■ ( x JVb, n\^-Na, T a ) , 

These one-particle amplitudes fulfill the semi-group property 

(x c ,T c |x a ,T Q ) = J d 3 X b (x c ,T c |x b ,T b ) (x b ,T b |x Q ,T Q ) 

and are invariant with respect to translations in the imaginary time: 

(x b ,T b +T |x a ,T a +T ) = (x 6 ,T b |x a ,T Q ) . 



(3) 



(4) 



(5) 
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We now impose the indistinguishability of particles. In D > 2 space dimensions there are only two kinds of indistin- 
guishable particles: bosons with a completely symmetric and fermions with a completely antisymmetric TV-particle 
wave function. According to Pauli's spin-statistic-theorem they are associated with integer and half-integer spins, 
respectively. The corresponding modifications of the path integral ([TJ) are straight-forward. For bosons, we have to 
sum over all permuted final positions 'X-phAb, where P(v) denotes one of the N\ possible permutations of the index v. 
For fermions, there is an extra minus sign for odd permutations. Restricting ourselves in this paper to the bosonic 
case, we must calculate the imaginary-time evolution amplitude 

(x lfc , . . . , X Nb , T h |x la , . . . , X Na ,T a ) B — — ^(Xp^jj,, . . . , Xp(jv)b, T b |x lQ , . . . , XjVa, T a ) . (6) 



III. CYCLE EXPANSION 



The canonical partition function of N non-interacting bosons is given by the trace 



7 b 



09) = J d 3 Xl ■ ■ • J d 3 X N (x 1 ,...,X N ,hp\Xi,...,X N ,0) B , (7) 

where (3 = (fcpT) -1 represents the reciprocal temperature. The calculation of this basic statistic quantity is non- 
trivial due to the indistinguishability of the particles. At first we consider the Bose gas with small particle numbers 
N . For N = 1, the partition function ([7]) reduces to the well-known expression 

Zxifi) = J d 3 x(x,hp\x,0) . (8) 

Here we have omitted the index B, as the distinguishability of particles is not an issue for a single particle. For N = 2, 
we obtain from ([7]) according to the indistinguishability © and the factorization property ([3]) 

Zi(f3) = ^ J d 3 x 1 d 3 x 2 {(x 1 ,h{3\x 1 ,0) (x 2 ,fr/3|x 2 ,0) + (x 2 ,n / 3|x 1 ,0)(x 1 ,n/3|x 2 ,0)}. (9) 

Using the translational invariance ([5]) together with the semi-group property ([U for the second term, this yields with 

m 

Zi(P) = \{zm+ Zi(2/3)}. (10) 
In a similar way, we obtain the result for N = 3: 

Zf(/3) = \{zl((i) + 3 Zi(/3) Zi(2p) + 2^(3/3)} . (11) 
For an arbitrary number N of bosons we must calculate according to and © the canonical partition function ([7]) 

Z M = ^J2j d3xi --- / ^(xpd,,^^) ••• (x P(JV) , ^1x^,0). (12) 

Due to the indistinguishability of the particles there are single-particle amplitudes contributing to the partition 
function (|12[) . whose final points coincide with the initial points of other particles. These amplitudes can be combined 
to many-particle formations, which can be represented by paths winding around a cylinder of circumference h(3. A 
typical situation for three particles is shown in Figure [TJ This represents a single closed cycle with three windings 
around the cylinder. An arbitrary single cycle with n windings is given by the multiple integral 

Z c n (P) = J d 3 x x - J d 3 s„(x 1 ,n/3|x„,0)(x n ,n/3|x„_ 1 ,0) ••• (x 3 , ty?|x 2 , 0) (x 2 , h0\x u 0) . (13) 

Due to the translational invariance ([5|) and the semi-group property (j4|) of the 1-particle imaginary-time evolution 
amplitude we see that 

Z%jS)=Z 1 {nP). (14) 
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FIG. 1: Example of paths contributing to (|13[) for a cycle of length n = 3: a) final points coincide with initial points of the 
other particles, b) and c) show the same situation in an extended zone scheme and wrapped upon a cylinder, respectively. 
Figure after Ref. )\^. 



Thus, a contribution of a closed cycle of length n corresponds to a simple one-particle partition function but with a 
temperature scaled down by the same factor n. With this we return to the general expression (jTSJ) for the canonical 
partition function. The full TV-particle partition function decomposes into mutually disconnected groups of closed 
paths winding around a cylinder of circumference hf3. Due to the indistinguishability of particles this decomposition 
is nontrivial since the cycles may have different winding numbers. 

In order to illustrate the situation, we consider in detail the possible path combinations for N — 2 and N = 3 
particles. Figure [2] shows the respective permutations in different notations. In the cycle notation the indices of 
particles contributing to the closed paths are stated within parentheses. This cycle notation is mathematically useful 
but is not adopted to the indistinguishability of particles as it relies explicitly on the numbering of the particles. 
For instance, the second, third, and forth permutation in the second row of Figure [5] are regarded as different 
permutations, although they are identical from the physical point of view. This motivates to introduce a different 
notation for the permutations which respects the indistinguishability of the particles. It is the cycle-number notation 
which characterizes each permutation by the 7V-tupel (Ci, . . . , Cn) of numbers C n of cycles with the length n. In this 
notation the second, third, and forth permutations in the second row of Figure [5] have the same 3-tupel (1, 1, 0). 

Using the result (fT4")) for the multiple integral (fT3"|) we arrive at the cycle decomposition of the canonical partition 
function (fT2"]): 

SC„n=JV 

P n=l 

Note that the numbers of cycles C n of length n, which occur for a certain permutation P, are restricted by the 
condition T, n C n n = N. The sum in (|15[) can be rearranged as follows. Instead of summing over permutations and 
decomposing them into cycles, we may sum directly over all different cycle numbers (Ci,...,C/v) which respect 
the equation £ n C„n = N. By doing so, we have to take into account that a particular configuration (Ci, . . . , Cat) 
can occur with a certain multiplicity M(C±, . . . , Cn). For instance, in Figure [2] we can read off M(3, 0,0) = 1, 
M(l, 1,0) = 3, and M(0,0, 1) = 2 for three particles. Thus, we obtain 

£C„n=iV oo 

Z n^) =m E M ^ ■ ■ ■ > C ^ II i z ^)f n . (16) 

C U ...,C N n=l 

where we have formally extended the product to all integers by setting C n = for any n > N due to the condition 
^nC n n — N. Now we derive a general formula for the multiplicities M(C\, . . . , Cn). In total the cycle numbers 
(Ci, . . . , Cn) can be arranged in N\ different combinations. However, not all of those represent different permutations 
as we shall demonstrate by the following examples: 



• The cycle numbers (Ci = 0, C 2 = 0, C 3 = 1) allow for 6 combinations (123), (312), (231), (132), (213), (321), 



5 



(1)(2) 
1 1 

2 2 

(2,0) 



(12) 
1 ^ 



2 

(0,1) 



(1)(2)(3) 
1 1 

2 2 

3 3 

(3,0,0) 



(12)(3) 
1 \ / 1 



(1,1,0) 



(1)(23) 
1 1 



(1,1,0) 



(13)(2) 

1 x / 1 



(1.1.0) 



(132) 

1 x / 1 



(0,0.1) 



(123) 



3 / \ 3 
(0,0,1) 



FIG. 2: Permutations for N = 2 (first row) and AT = 3 (second row) particles. The indices above each permutation stand for 
the cycle notation, whereas the respective iV-tupel (Ci, . . . , Cjv) below each permutation denotes the cycle number notation. 



where the former and latter three combinations are different representations of only two permutations 



3 12/' ' V 2 3 1 , 

• The cycle numbers (Ci = 3, C 2 = 0, C 3 = 0) lead to 6 combinations (1)(2)(3), (1)(3)(2), (2)(1)(3), (2)(3)(1), 
(3)(1)(2), (3)(2)(1) which are all representations of the same permutation 

These observations are generalized as follows: 

• A cyclic permutation within a n-cycle does not lead to a new permutation, so we have n irrelevant combinations. 

• Exchanging cycles of the same length does not lead to a new permutation, so we have C n \ irrelevant combinations. 

Summarizing these combinatorial considerations, we obtain for the multiplicity the Cauchy formula of the permutation 
group 



2l|,l22|: 



Af(C 1 ,..,C,) = n -^. (19) 



After this combinatorial excursion we return to the canonical partition function (I16[) . Using (|19p . we find our second 
main result in form of a cycle representation of the canonical partition function [lO} : 



£C„n=JV oo 



*%w= e n?fi 



c 

Ci,...,C w n=l " 



Zi(n(3) 



c„ 



(20) 



Note that specializing l[20|) for N — 2 and N = 3 particles reproduces, indeed, the earlier findings ([10]) and (fTTj) . 

Although the cycle representation allows us, in principle, to determine the canonical partition function Z^(f3), 
this approach is not efficient for a large particle number N as the computational effort grows factorially with N. 
Furthermore, it is rather cumbersome to figure out all cycle numbers C„ which respect the condition Y, n C n n = N. 
Therefore, we need to find a more efficient algorithm to determine canonical partition functions. This is done by using 
the fact that the grand-canonical partition function is the generating function of all canonical partition functions []~9l ] . 

IV. GENERATING FUNCTION AND RECURSION RELATION 

Introducing the fugacity z = e 13 ^ with the chemical potential /i, the grand-canonical partition function reads 

00 

ZS c (!3,z)=J2Zn(P)z N - (21) 

N=0 
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This quantity serves as the generating function of canonical partition functions Z^(/3), which are simply the Taylor 
expansion coefficients 



1 d N Z§ c ((3,z) 



dz 



N 



(22) 



z=0 



In order to evaluate (|2T|) , we use the cycle representation (|20p and insert the condition N — YjC n n into the exponent 
of the fugacity z. Thus, we obtain at first 



oo T,nC n =N oo 



N=0 d,...,C N n=l 



Cm 



Z 1 (nj3)-A 



Ci=0C 2 =0 n=l 



Cn 



(23) 



In the second expression the condition T,nC n = N no longer occurs. Due to this simplification we can rewrite the 
multiple sum over the cycle numbers C n as a product over the individual sums and obtain 



oo oc 



-zfcOM=n e 



71=1 C„=0 



1 

CJ. 



Zi{n(i)z n 



(24) 



Performing the C^-summation yields the result 



Zf c (/?,z)=expi£ 



Zi(n(3)z 



(25) 



The partition function of the vacuum, where no particle is present (N = 0), is obviously found from (|22p with (|25|) : 

Z${p) = l. (26) 
For N > 1, we observe that the grand-canonical partition function ([25]) obeys the differential equation 

(27) 



dzBc d f z) =zB c (p, Z )jrz 1 mz k -K 



k=l 



This is a self-consistent equation for the generating function which we solve recursively. Forming the iVth derivative 
of (|2T1) via the Leibniz rule of differentiation: 



d N Zg c ((3,z) 



dz 



N 



z=0 



N 

E 



< (n- 1)! (JV-n)! 



(9.3 



z=0 u= 



J2z 1 (kp) 



Qn—l^k—l 



k=l 



dz r ' 



z=0 



and taking into account the identity 



Qn—l^k—l 



dz n ~ 



(n - l)\8 k ,n , 



(28) 



(29) 



z=0 



we obtain from (|2"2")) the desired recursion relation for the canonical partition functions Z^{(3) [EEl, 3 IS, 20 1: 

1 N 

Z§(/3) = ^ E Z ^ Z N-n(P) ■ (30) 

n=l 

Inserting the starting value (f2"6")l . we immediately reproduce for N — 2 and TV = 3 the earlier results (fit))) and 

For large particle numbers it can be shown that the computational effort grows only algebraically with N. Thus, 

iterating Eq. ([50)1 can be done very efficiently, as is further worked out in Ref. [l8| . 

Once we have obtained the ./V-particle partition function, we can proceed to calculate all thermodynamic quantities 
of interest. For instance, in Section IVlII we shall determine the specific heat capacity according to 



C*(T) = k B T-^[T In Z»] 



(31) 
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V. SPECTRAL DECOMPOSITION 

In order to apply the recursion relation (|30j) to specific systems we must calculate the 1-particle partition function 
Zi((3) of Eq. ©. For this we solve the time- independent Schrodinger equation and use the energy eigenvalues E^ to 
calculate Zi((3) as a spectral sum 

Zi 08) (32) 

k 

where k stands for the quantum numbers. The n-fold cycle (|14[) 

Z c n ((3)=J2^ (nEk) (33) 



may be interpreted physically as a partition function (|32| of a single quasi-particle consisting of a correlated union of 
n fundamental particles with the total energy nE^. 

The spectral representation (|32|) leads via (|25|) to the grand-canonical partition function 



Z* C (J3, n) = exp J - £ In [l - e^"*)] I . (34) 



k 

From it we may derive the average particle number in a grand-canonical ensemble according to 

{N) = \ | Pa*&(M] = E ^) _ i - (35) 
which is the well-known Bose-Einstein distribution. 



VI. GROUND-STATE OCCUPANCY 

The most important physical quantity in Bose-Einstein condensation is the condensate fraction iVo which is also 
called the ground-state occupancy. In the usu al g rand-canonical treatment this quantity is found by solving numer- 



ically the Bose-Einstein distribution l|35p [19l . |20( . To this end one calculates the chemical potential fi for a given 
temperature T at a fixed average particle number (N). Given fi(T) one finds the average number of particles in 
ground-state fraction from the equation (No) — le^^ 0-1 ^ — l] , where Eq is the ground-state energy. 

In the present paper we shall treat this problem canonically at a fixed particle number N, without fixing a chemical 
potential. To this end we calculate directly the probability that one of the bosons resides in its ground state [15|, ll6| . 
This probability multiplied by N will be considered to be the finite- N version of the condensate fraction ./Vq. In order 
to determine this quantity, we isolate in Section IVI Al the ground-state contribution to the cycle representation (|20p • 
Subsequently, in Section [VI Bi we derive a more efficient equation for the canonical ground-state occupancy. 

A. Cycle Representation 

Using the spectral decomposition (|32|) . we can rewrite Z\(n(3) according to 

Zi(n/3) = 7«G9)+£«09). (36) 
The first term denotes the contribution of the ground state 

7 „(/3)=e-"^ G , (37) 
where Eq denotes the ground-state energy. The second term in (f3"o| stems from all excited states 

£„(/?)= ]T' e-"^", (38) 



8 



where the prime indicates the omission of the ground state. Using the binomial formula, the canonical partition 
function (|2"U)) is then decomposed as follows: 



T.nC n =N Ci Cjv 

zm= E E - E - 

Ci,...Cjy mi— mjv— 



E nm n 



N 



e. 



C„-m„ 



09) 



li - m„! (C„ - m„) 

n— 1 v 7 



lrc.C„ 



(39) 



where we have used the factorization property of the ground-state contribution (|37p . which reads j n {P) — 7"(/9)- 

Now we are able to determine the weight of the particles residing in the ground state. Consider first the cases of 
N = 2 and N — 3 particles, where (j39|) yields 



z 2 B = 7i°^ 1 2 +6]+7 1 1 £i+7i\ 



-1 



+ & +7f6+7i 



(40) 
(41) 

The number of 



This can, of course, be directly obtained from (fTU)) and (TTTj) by inserting the decomposition 
particles in the ground state coincides always with the power of 71. Thus, the ground-state weights and 
associated with and Z$ are 



= 0x 7 «- 



1 



1 



+ 1 x 7} a + 2 x 7l 2 



1 



1 1 



W 3 =0x7^ ^-^ + -^^ + -^ 3 j +lx 7 f- [£ + &J +2x7i^i +3x 7 ^ 
The general formula for the weight of the particles in the ground state is therefore given by 

SnC„=« d C N / N \ N tC„^m„ 



wgw = E E - E E nm n [7i(/3)] Enm - n 



Ci,...Cn mi— mjv— \n— 1 



J- 1 ; m n \ {C n -m n )\n' 



(42) 
(43) 

(44) 



B. Recursion Relation 



The evaluation of the cycle representation (|44|) of the ground-state weight is quite cumbersome for large particle 
numbers N. Thus, we aim at deriving a more efficient relation for determining W^{(3). Therefore, we observe that 
the canonical partition function |39|) and the ground-state weight (|44|) are related via 

d 



WW) =7i 



-zm- 



dim 

It is useful to define a generating function for the ground-state weight as the sum 



(45) 



from which we conclude with Eqs. (f2"Tj) and (j4"S"|) 



>V| c (/3^)=7i(/9) 



5 



5 7 i 09) 

Taking into account the decomposition (|36p we see that according to 

oo 

W| c (/3, z) = Z| c (/3, z) E 7n(/3) • 
n=i 



(46) 



(47) 



(48) 



The canonical ground-state weights can then be found according to Eq. (|46[) as the Taylor coefficients following from 



1 d N W§ c (J3,z) 



Nl 



dz N 



z=0 h 



N 1 

TuW. 



d N - l ZS c ((3) 



^ll(N-l)\ dz N - 1 



E 7n09) 

2 = „ = 1 



(49) 



z=0 
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Using Eqs. ([22)) and (|29| . we obtain the canonical ground-state weights as 

N 

W^((3) = Y / ln((3)Z§_ n ((3), (50) 

n=l 

which represents the result of Ref. [l5| . It states that the weight W^(f3) of the particles residing in the ground state 
of an iV-particle Bose gas is immediately known, once the canonical partition functions Z^((3) are calculated for all 
n<N. A proper normalization of the statistical ground-state weight W^((3) allows us to determine the probability 
wfj{f3) to find a particle in the ground state: 



This quantity multiplied by N represents the condensed fraction Nq of an TV-particle bosonic ensemble. 



VII. SPECIAL POTENTIALS 

Now we specify the potential energy V(x) and calculate the heat capacity and the ground-state occupancy within 
canonical ensembles. In order to see what can go wrong, we first follow Feynman's textbook on Statistical Mechanics 
[l(J and describe in Section [VII Al the homogeneous Bose gas, which leads to unphysical results. This is corrected in 
Section lYlIB I bv a modified approach for the homogeneous gas where the special role of the ground state is properly 
taken into account. This approach corresponds to the semiclassical treatment of a Bose gas in a box- like trapping 
potential, which is then investigated quantum mechanically exact in Section IVII CI 

Although, the quantities of interest may be calculated in principle analytically from Eqs. ([50)1 and ([ST))) , the calcu- 
lations in this section are performed numerically with a C-program in order to reach large particle numbers. More 



details concerning analogous numerical implementation techniques are presented in Ref. 19 1. 



A. Homogeneous Bose Gas 

Consider a gas of non-interacting bosons without any trapping potential, V(x) = 0, i.e. the homogeneous Bose 
gas. This theoretical model represents a paradigm studied in the seminal work by Bose and Einstein [23]. Also the 
original work of Feynman [Io| within the canonical ensemble relies on treating the homogeneous Bose gas. But this 
historical important model has problems, which manifest themselves in a divergent ground-state occupancy and in a 
violation of the third law of thermodynamics. To exhibit this explicitly, we insert the free-particle energies 

£ k = — k 2 (52) 

with continuous wave vectors k into the one-particle partition function (|32[) and perform the sum over all energy levels 
in form of a continuous integral 



with V being the volume of the system. The result is 



Z iM) = J^> ( 54 ) 



where A = y/ '2ir (3h 2 / 'M denotes the thermodynamic de Broglie wave length. According to ([52]) the ground-state energy 
vanishes for the homogeneous case, thus leading to the ground-state contribution ([57)) of the partition function 



7i(n/3) = l, 



(55) 
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FIG. 3: a) Heat capacity per particle and b) ground-state occupancy versus reduced temperature in canonical ensembles of 
N =1, 2, 10, 100, 1000 particles (dashed lines) compared with the grand-canonical result in thermodynamic limit (solid line). 
The particle density is chosen to be constant. 



which is independent of temperature. Now we use a dimensionless temperature parameter 



_ y2/3 W rp 

2ttH 2 



(56) 



and represent our results for one and two particles. From recursion relation ([50)1 we obtain the canonical partition 
functions 



Z?{r) 



3/2 



r 3/2 
25/2 



T 

y 



(57) 



which also follow from (fTU)) . (fTTj) with (|51)) . From (|57[) we calculate the heat capacity according to the thermodynamic 
equation (|31|) . Using again the dimensionless temperature parameter (|56[) yields for one and two particles 



C?{t)= -k B , 



C 2 s (r) = 3fci 



1 + 9y/2T 3/2 + 16t 3 

2 + 8V2t 3 / 2 + 16 t 3 



(58) 



From both results we read off that they obey in the high-temperature limit the classical Dulong-Petit law C B 



N 



3-/Vfce/2- However, in the small-temperature regime one obtains the general result C B — > 3fce/2. This finding 
violates the third law of thermodynamics, according to which the heat capacity has to vanish at zero temperature. 
Both results (|55|) are plotted in Figure [3] a) versus the reduced temperature t = T/T c , where 



T c = 



2irh 2 
k B M 



N 



VC(3/2) 



2/3 



(59) 



represents the critical temperature of an ideal homogeneous Bose gas. The reduced temperature t is related to t of 
Eq. $5U) by 



t = 



C(3/2) 



N 



2/3 



(60) 



In order to obtain results for larger particle numbers, we evaluate the recursion relation (|3"D|) . Here the variable (|56"]) 
is more useful than the reduced temperature (|60[) since the latter depends explicitly on the particle number. The 
recursion (|30p involves results for different particle numbers, so that the employment of the reduced temperature t 
would require a permanent rcscaling. Plotting our results for different particle numbers in Figure [3] a) then only 
requires a rescaling of the original variable r to the reduced temperature t according to (|6H)) . This plot shows the 
behavior of the heat capacity per particle for different total particle numbers. Here, the well-known grand-canonical 
result is reached in the thermodynamic limit. Furthermore, the heat capacity per particle is suppressed at zero 
temperature in the thermodynamic limit, but, as mentioned before, the heat capacity itself remains non- vanishing for 
t = 0. 
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A related but much more dramatic problem arises for the homogeneous Bose gas in finite systems when we consider 
the probability to find a particle in the ground state. To illustrate this we calculate first from (1ST))) the ground-state 
weights for one and two particles using ([5l 



W*(t) = 1, W 2 B (r) = l + r 3 / 2 , (61) 

which also follow from (|4"2")) and l)4"3")) . The ground-state probability (|5"T[) then yields with (|57[) and (f6"Tj) 

td . . 1 R/ . 2V2 1 + t 3 / 2 . , 

^=^> ^^ l + 2V 2,3/ 2 - ^ 

Here we see that the probability to find a particle in the ground state decays for increasing temperature but diverges 
for low temperatures. This behavior is depicted in Figure [3] b) where the ground-state occupancy is also shown for 
larger particle numbers. One can also see that our canonical results do not match the well-known grand-canonical 
shapes in the thermodynamic limit. Moreover, it diverges for all finite systems at zero temperature. Not only that 
this divergency is physically meaningless, but by definition any probabilities have to be bounded by unity. In the 
present description of the homogeneous Bose gas in the canonical ensemble this is obviously not the case. As it will 
be seen in the next subsection, this problem appears due to an improper treatment of the semiclassic limit and the 
results can be corrected by taking care of the special role of the ground state. 



B. Homogeneous Bose Gas with Ground State 

Inspecting the results in the previous subsection more carefully, we recognize already at the beginning a contradiction 
between the one- particle partition function and its ground-state contribution (j5"5)) . Namely, whereas ([M]) vanishes 
for low temperatures, the ground-state partition remains finite. This leads to a negative contribution for all excited 
states. However, one would reasonably expect that the latter has to be always positive and vanishes only for zero 
temperature. According to this expectation we suppose that the partition function has to be a sum of the ground-state 
contribution (|55[) and the contribution of excited states in the form given by ([51 



Zi(nj9) = 1+ . (63) 

Such a decomposition is equivalent to replace condition ()53[) by 

£/(k)^/(0) + F| /(k) (64) 

for any function /(k). Comparing (|64|) with relation (|53[) we see that the ground state in the latter didn't play any 
special role. There we integrated over all states including the ground state, so that the latter was only a point of zero 
measure, as basically any other state. In our present approach, however, we separate the ground state explicitly from 
the remaining excited states. In order to demonstrate how well this procedure works, we calculate the heat capacity 
(f3Tj) of one particle from ([63]) . yielding for D = 3 spatial dimensions 

B 3k B 5r 3 / 2 + 2r 3 

^ (T) = — 2 + 4r3/ 2 + 2T 3 • (65) 

This leads to 3fcs/2 in the high-temperature limit according to the Dulong-Petit law. For small temperatures the 
heat capacity ((65)) vanishes satisfying also the third law of thermodynamics. Both results (|58|) and (|65|) for the heat 
capacity of one particle are compared in Figure |4] a). 

The second problem of the homogeneous Bose gas in the last subsection concerns the probability to find a particle in 
the ground state which diverges for vanishing temperature. Taking into account the extra ground state, we recognize 
that (|55p represents its contribution to the partition function (|63[) . Thus, the probability for a particle to be in the 
ground state seems to be normalized. For one particle we use Eq. (|50|) to obtain W^(t) — 71 (r), yielding with the 
definition (|5Tj) the ground-state probability 

<M = ttW • (66) 
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FIG. 5: a) Heat capacity per particle and b) ground-state occupancy versus reduced temperature calculated with an extra 
ground state according to ()63|) . Results are represented for canonical ensembles of N = 1, 2, 10, 100, 1000 particles (dashed and 
dashed-dotted curves) and for the grand-canonical ensemble in thermodynamic limit (solid curve) for a fixed particle density. 



It vanishes for high temperature (t 3> 1) and has its maximal value of one for zero temperature, as one may expect on 
general physical grounds. We compare this result with the corresponding one of the previous subsection in Figure [4] 
b). For small temperatures the corrected curve does not diverge and shows the expected behavior. One understands 
immediately, why this procedure works so well for finite systems. Assuming the particle density to be equal for all 
systems, the volume V has to be small for small particle numbers. The additional unity in (|63p plays for them a more 
important role than for larger ensembles. Moreover, this unity is practically negligible for high enough temperatures 
and carries a larger weight for lower temperatures. As the contributions of different lower particle numbers are mixed 
in the quantities of interest (|30f and (|50|) . this influences the situation decisively. 

Now we turn to many-body systems and calculate their heat capacity according to the recursion relation (|30|) 
with definition (|3Tj) and the ground-state probability according to (|50[) with (|5Tj) . The corresponding results are 
represented in Figure O All curves have the correct behavior for high temperatures as well as for the low-temperature 
region. Moreover, in the thermodynamic limit, where the particle number goes to infinity, both results display the 
grand-canonical shape. This means that, although all curves are smooth for finite systems, the behavior of the heat 
capacity and the ground-state occupancy change more pronounced at the critical point by increasing the number of 
particles. 

Thus, we can summarize that a careful treatment of the ground state in this subsection improves the description 
of the homogeneous Bose gas in a fundamental way. Nevertheless, also this attempt represents only a physical 
approximation to a real finite system. To this end we will show in the next subsection that the exact quantum 
mechanical results deviate quantitatively in a significant way from the semiclassical treatment presented here. 
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C. Bose Gas in Box Potential 



In this section we work out the exact quantum mechanical treatment of a Bose gas in a square-box potential. For 
convenience, we choose the size of our box to be equal to L in all spatial directions, so that the trapping potential 
reads 



V(x) 



, 0<Xj<L for all j = 1,2,3, 
oo , elsewhere . 



The one-particle Schrodinger equation is then solved exactly by the energy eigenvalues 

7T 2 fl 2 . 



2 ML 2 



with kj = 1,2,3,4,... (j = 1,2,3) 



(67) 



(68) 



As one can deduce from the normalization condition of the eigenstates, the value k = (ki, &2, fe) is not a solution, if 
at least one kj vanishes. Thus, the ground-state energy is given by Eq = 3tt 2 H 2 / (2ML 2 ). From this we obtain the 
ground-state contribution to the partition function (|37l) which reads with the dimensionless temperature parameter 
(156]) as 



7i(n/3) 



_ e -3n7r/(4r) 



The one-particle partition function then reads 



E 



-rnvm /(4t) 



(69) 



(70) 



The corresponding quantities in the Subsections I VII Al and lVII Bl are different semiclassical approximations of (|69|) and 
(I70p . which are justified for an infinitely large box. Such a semiclassical approximation leads, at first, to the ground- 
state contribution (I55|) and, secondly, to a vanishing level spacing which gives rise to the continuously distributed 
energy levels (|52|) . Therefore, the sum in ([70]) has to be replaced by an integral according to (|53[) or (|64[) . respectively. 
The one-particle partition function is then given by a Gaussian integral and can be simply calculated leading to (|54p 
or (|63[) . Strictly speaking, this semiclassical description of a finite system is inadequate insofar as it is impossible to 
have a finite number of particles in an infinite box, unless in the trivial case of zero density. 

The partition function (|70[) in a box cannot be calculated exactly. Thus, we can evaluate the sum only approxima- 
tively up to a certain order of the summation index m. However, this procedure does not work well for all temperatures 
since the exponential functions e - n7T /( 4l ~) i n (|70| are required to be small, and this is only true for small temperatures 
(r <C 1). In order to perform the sum in (I70j) also for high temperatures, we need its dual expression, which follows 
from the Poisson summation formula 11 | 



E /( m ) = E / dzf(z)e 2 ^ 



(71) 



q— — oo 



With this, the one-particle partition function ([70)) is converted into 



Zi{n[3) 



V 9=1 



(72) 



Of course, this function cannot be calculated exactly either, so we have to perform the summation also here only 
up to a certain order. Whereas the approximation (|70|) is applicable for small temperatures, its dual ([72]) can be 
evaluated for high temperatures. How can we now use both representations (|70l) and (|72[) in order to evaluate the 
one-particle partition function in the whole temperature range? In principle, both representations are equivalent, but 
only when both series are evaluated exactly. This exact expression is approximated well by summing few terms in 
(|70p only below a certain temperature and in (|72|) only above a certain temperature. There exists an appropriate 
temperature region where both approximations coincide. In Figure[H]we see that for the interval 0.1 < r/n < 0.5 even 
the lowest term in (|70l) as well as (|72[) are sufficient. By taking more terms into account, this matching region grows. 
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FIG. 7: a) Heat capacity per particle and b) ground-state occupancy versus reduced temperature in exact square-box poten- 
tial for ensembles of TV = 1,2, 10, 100, 1000, 10000 bosons (dashed and thin curves) and in the grand-canonical ensemble in 
thermodynamic limit (thick solid curve). All results are calculated for the same particle density. 



For practical calculations, however, the summation of many terms costs much effort, so it is more efficient to precisely 
define an intermediate temperature, which joins the low-temperature region with the high-temperature region best. 

In order to calculate the results for many particles we use the recursion relation (|30[) for the partition functions and 
({50")) for the ground-state weights. The heat capacity is then calculated from (|31[) and the ground-state occupancy 
from iJSTJ). Results for N = 1,2,10,100,1000 and 10000 bosons are depicted in Figure [7] All results for the heat 
capacity exhibit the Dulong-Petit behavior in the high-temperature limit and fulfill the third law of thermodynamics 
for T — > 0. The ground-state occupancy is 100% for zero temperature and vanishes for high temperatures. One specific 
feature of an exact quantum mechanical treatment is worth emphasizing. Due to the non-vanishing gap between the 
ground state and the first excited state, the slope of the heat capacity is exponential for small temperatures. This 
is especially visible for small particle numbers TV. For higher TV, the exponential region decreases and is followed by 
the usual power-law increase ~ t 3 / 2 . In the thermodynamic limit the exponential region disappears and the power 
behavior takes over. 

As expected, the canonical results converge in the thermodynamic limit towards the grand-canonical one for an 
infinite average particle number. However, this convergence is much slower in the quantum-mechanical exact treatment 
than in the semiclassical approximation of the previous subsection. This manifests itself, for instance, in an apparent 
discrepancy in the location and extent of the heat capacity maximum even for TV — 10000 particles. The same 
observation holds for the ground-state occupancy. Although an increase of the particle number makes the curves 
sharper in the crossover region, the location of the latter is clearly shifted towards higher temperatures. This gives 
rise to the so-called finite-size effect which will be elaborated in the subsequent section. 
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VIII. FINITE-SIZE EFFECT 



Let us now investigate in detail the crossover temperature, which can be regarded as a generalized critical tem- 
perature T c . At first, we suggest to define this as the temperature where the ground-state occupancy possesses the 
largest curvature. This definition is motivated by the fact that in the thermodynamic limit this quantity has an 
infinite curvature at the critical point. By analogy, in a finite ensemble the point with a significantly sharp rise in the 
curvature may play the role of the transition temperature. These points are plotted with dots in Figure [8] for different 
particle numbers. The finite-size effects are quite large. In a box confinement, they amount to 30% for TV = 300 and 
to about 5% for N = 300000 bosons. The effects here are about 3-4 times larger than for the Bose gas confined in a 
harmonic trap 0|. Furthermore, it is remarkable that the finite-size shifts of T c point upwards which is opposite to 
the harmonic case. This fact can be explained by a stability analysis, which shows that in a 3D-box potential only 
finite systems are stable ■ Particle fluctuations grow with the size of the box, resulting in a loss of stability, and 
this is why bigger systems condense at lower temperatures. 

Let us compare our canonical results for the critical temperature with analytic approximations. The first estimate 
of the shift of the critical temperature was given in Ref. [25| by 



AT C = 1 

T c C 2/3 (3/2)iW3 n 



4 AT 2 / 3 



3ttC 2 / 3 (3/2) 



(73) 



This result is plotted in Figure[8jas a dashed curve. It represents the grand-canonical result for semiclassical corrections 
caused by large but finite particle numbers N, which takes into account the non-vanishing energy eigenvalues kj of 
(f68|) in a box potential. This correction was calculated in Ref. by replacing the sum over all energy eigenstates 
in (|35| by a momentum integration and subtracting all states lying at coordinate planes in the momentum space. As 
one can see in Fig. [8j this analytic result reproduces well the direction of the critical temperature shift but drastically 
underestimates its magnitude. Further corrections could be obtained by avoiding double-countings at cross sections 
of coordinate planes as well as by taking into account the finite gap between the ground state and the first excited 
states. 

Here we prefer a more controllable procedure for calculating the finite-size corrections which goes as follows. We 
start with the cycle number representation of the grand-canonical average of the particle number according to Eqs. 
and (1331): 



(N) = Zi(n(3)z n = (N ) + (N T ) ■ (74) 

71=1 

It contains the number of particles in the ground-state condensate 

W = i;e-**V= z -i e ,E G _ 1 (75) 
and the number of thermal particles in excited states 

oo 

(N T ) = i Z i M - ^ n0EG \ z n ■ (76) 

n=l 

Furthermore, we rewrite the one-particle partition function ([70]) as 

Z 1 (n/3) = e~ 3 "^ 4r ) [l + e - 3 ™ 7r/(4T V(r/n)j 3 (77) 
with dimensionless temperature (|56|) and an excitation contribution 

oo 

CT (r) = e-^-^l^ = ^ - 3/2 + vr/V^ + 0(l/r) , (78) 



m=2 



where the latter equation is the high-temperature approximation similar to the duality relation (|72[) for the whole 
partition function. Since we aim at finding the location of the critical temperature, we demand for the fugacity z the 



16 



condition ze i3Eg 



ze 



-3tt/(4t) 



1 but at the same time (No) <C (N). Both conditions are exactly fulfilled in the 



thermodynamic limit and are reasonable approximations for large but finite systems. From Eqs. (|74p - (|78p we obtain 
at the critical point t = r c approximately 



oo 

W = £ [3e~ 3n ^^a(Tjn) + 3e- 6n ^^<j 2 (T c /n) + e- 9n */^<r 3 (T c /n) 
=1 



(79) 



According to Eq. ([78]) the sums in ([79]) can be expanded in terms of the polylogarithmic Bose functions C, v (e x ) 



(N)=r 3 J\ V2 {e-^'^) + 



3r e 



+ ■ 



3r, 



1/2 



!Ci (e^/^)) -3Ci(e- 9V(4rc) ) 

'4Ci/a (e- 8 */^)) - 12& /2 (e- 6 ^)) + (4tt + 9) Cl/2 (e" 9 ^)) 
For large values of r c , the polylogarithmic functions are evaluated with the help of the Robinson formula 26 1 

oo 

T(l - v) x"- 1 + (- x ) k C( v -k)/k\ for 1 v < 2 , 

oo 

-\nx+^(-x) k ((l-k)/k\ for v = l, 



(80) 



(81) 



k=l 



where T(z) denotes the Gamma function and £(z) = X^^Li I/" 2 the Riemann zeta function. Inserting this into 
and identifying (N) with the actual particle number N , we arrive at 



N 



r 3/2 



C(3/2) 



In 



Co 



2t c 



3(tt+1) 



-y 2 C(l/2) + O(r c ) 



Using the reduced temperature t = T/T c instead of r according to ([6"0")) , this can also be rewritten as 



.3/2 



2 iV 1 /3^2/3(3/ 2 ) 



111 



2N 2 ' 3 t c 
ttC 2 / 3 (3/2) C~ 3 



3(vr+l) ic /2 C(l/2) 



JV2/3 ^1/3(3/2) 



0(1/N) . 



(82) 



(83) 



Taking into account only the first term on the right-hand side yields the semiclassical result t c = 1. The first correction 
involves a coefficient C3 which unfortunately cannot be calculated directly. The origin of this problem is the divergency 

of the polylogarithmic function Cv( z ) f° r z ~ 1 an d any v < 1, which is explicitly given by the first term in (|81[) . So 

1/2 

one can see that the last three terms in Eq. (|80j) contribute not only to the order of t c but also to the order r c . 
The same is valid for all other terms which we omitted in ([50)) . The leading finite-size correction of r c plays here a 
special role insofar as it consists of infinitely many terms and thus its factor C3 cannot be calculated from the first 
few terms. However, we can perform a direct numerical evaluation, yielding the value C3 ~ 0.9574 which also follows 
from a semi-analytic calculation 11] . In contrast to that, all higher corrections in ([55)) are analytically calculable. 

In fact, such a problem is not typical for the box potential only. In a harmonic trap, one encounters an analogous 
situation for the subleading finite-size correction instead of the leading one for the present box case. Here, we should 
mention, that it is just this non-analyticity of the leading finite-size effect, which is responsible for the failure of 
the approximation in Ref. [25]], see Eq. ([75)) here. Their perturbative approach yielded the factor 3/2 instead of the 
numerical correct C3 -value. 

By extracting the A^-dependence of the reduced critical temperature of Eq. (|83p. we now obtain a grand-canonical 
semiclassical result which is represented by the solid line in Figure [5J As one can see, it overestimates canonical 
results considerably. In fact, although many physical properties within canonical and grand-canonical ensemble are 
equal in the thermodynamic limit, they can be quite different for finite systems. Let us now determine crossover 
temperatures for finite grand-canonical ensembles in the same way as before in the canonical approach. To this 
end we calculate the ground-state occupancy for a grand-canonical system. As already outlined at the beginning of 
Section IVI1 this is done by determining at first the fugacity z which corresponds to a given average particle number 
(TV) for a fixed temperature from ([74)) and inserting its value into ([75)) . The resulting quantity (Nq)/(N) is a smooth 
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FIG. 8: Shifts of critical temperature for N bosons in square-box potential with N = 300, 1000, 3000, 10000, 30000, 100000, 
300000 (from right to left) in canonical ensemble (dots). They differ considerably from the grand-canonical results (triangles) 
with the same boson numbers. The dashed curve corresponds to the semiclassical result (|73[l within grand-canonical ensemble 
from Ref. [25|]. The solid curve represents our result AT C /T C = t c — 1 according to Eq. f|S3[) . 

function of temperature similar to the canonical results in Fig. b) but somewhat shifted to the right. The crossover 
temperatures are then calculated by analogy with the canonical case as points of largest curvature and are indicated 
by solid triangles in Figure [5J They agree very well with the analytic approximation of (|83[) even for smaller systems 
of a few hundred particles, but differ strongly from canonical values. 

It is well-known that the canonical partition function can be obtained from grand-canonical partition function by 
means of the integral over the chemical potential fj, using the identity 

Z B n=P Se-^Z&OO. (84) 

In the leading saddle-point or steepest-descent approximation one obtains 

Z^ = e-^ N Z§ c (p)+ (85) 

where p, is the extremal point of the exponent — (3fiN + log Zq C (h) in (|84[) . The extrimity condition implies by means 
of Eq. ([34]) immediately the usual Bose- Einstein distribution ([35]) . For infinite N, this approximation becomes exact, 
leading to the equality of thermodynamical properties [l], \2l\ in canonical and grand-canonical ensembles. 

For finite N, one may want to calculate corrections from the neighborhood of the saddle point. The naively 
perturbation up to first order, howewer, yields a wrong sign of the shift for the critical temperature. A detailed 
analysis shows that a non-analyticity occurs at the saddle point near the critical point and infinitely many orders 
in the saddle point expansion have to be calculated and resummed to find the correct result. This fact has been 
discussed in literature, for instance in [l[ for the box potential and in [28] for a general case near the critical point, 
and is not yet solved. 

IX. CONCLUSION 

By taking into account the special role of the ground state of the system, we have obtained a consistent semiclassical 
description of an ideal Bose gas of a fixed finite number of particles confined in a box. We have thus improved 
Feynman's original canonical approach for treating the ideal homogeneous Bose gas. Finite-size effects have been 
shown to be significant, even for a moderate number of particles. Our canonical findings are compared with grand- 
canonical ones, for which we have found analytic expressions which improve results in the literature. An analogous 
analytic expression within the canonical ensemble could not be found, so we have to be content with numerical results, 
which were given in the present paper. 

Our canonical investigation may serve as a theoretical basis for future experiments in box-like traps manufactured 
with strongly inhomogeneous fields or laser beams. 
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